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Binary black-hole systems with spins aligned or anti-aligned to the orbital angular momentum, and which 
therefore do not exhibit precession effects, provide the natural ground to start detailed studies of the influence of 
strong-field spin effects on gravitational wave observations of coalescing binaries. Furthermore, such systems 
may be the preferred end-state of the inspiral of generic supermassive binary black-hole systems. In view of 
this, we have computed the inspiral and merger of a large set of binary systems of equal-mass black holes with 
spins parallel to the orbital angular momentum but otherwise arbitrary. Our attention is particularly focused on 
the gravitational-wave emission so as to quantify how much spin effects contribute to the signal-to-noise ratio, 
to the horizon distances, and to the relative event rates for the representative ranges in masses and detectors. As 
expected, the signal-to-noise ratio increases with the projection of the total black hole spin in the direction of 
the orbital momentum. We find that equal-spin binaries with maximum spin aligned with the orbital angular 
momentum are more than "three times as loud" as the corresponding binaries with anti-aligned spins, thus 
corresponding to event rates up to 30 times larger. We also consider the waveform mismatch between the 
different spinning configurations and find that, within our numerical accuracy, binaries with opposite spins 
Si = —S2 cannot be distinguished whereas binaries with spin Si — S2 have clearly distinct gravitational- 
wave emissions. Finally, we derive a simple expression for the energy radiated in gravitational waves and find 
that the binaries always have efficiencies E la .d/M > 3.6%, which can become as large as _E ra d/Af ~ 10% for 
maximally spinning binaries with spins aligned with the orbital angular momentum. These binaries are therefore 
among the most efficient sources of energy in the Universe. 

PACS numbers: 04.25.Dm, 04.30.Db, 95.30.Sf, 97.60.Lf 



I. INTRODUCTION 

It has been a long-standing goal of the field of numerical 
relativity to provide results for gravitational-wave data anal- 
ysis and thus enhance the capabilities of current and future 
gravitational wave detectors, in particular regarding the ob- 
servation of compact binary coalescence. With a series of 
breakthroughs in 2005 fl-IH], this long-term goal has suddenly 
become reality. However, much further work is required to ac- 
tually understand the practical implications of numerical so- 
lutions of the full Einstein equations for gravitational-wave 
data analysis. Indeed, first studies suggest that template banks 
that use numerical information can increase the reach of de- 
tectors a id the calibration of search pipelines 0-0], and 
improve the estimation of parameters, such as e.g., sky loca- 
tion or 



In this paper we use gravitational waveforms from 
numerical-relativity (NR) calculations of a number of se- 
quences of equal-mass spinning black-hole binaries whose 
spins are aligned (anti-aligned) with the orbital angular mo- 
mentum, and consider the detectability of these binaries for 
the ground-based gravitational wave-detectors as well as for 
the planned space-based LISA interferometer. 

Our interest in this type of binary stems from the fact that 
there are indications they represent preferred configurations 
in nature, at least if the black holes are supermassive. It has 
been shown, in fact, that when the binary is surrounded by a 
massive circumbinary disc, as the one expected by the merger 
of two galaxies, the dissipative dynamics of the matter pro- 
duces a torque with the effect of aligning the spins to the 



orbital angular momentum fllll . In addition, the merger of 
binaries with aligned spins yields recoil velocities which are 
sufficiently small (i.e., < 450 km/s OH-O) to prevent the fi- 
nal black hole from being expelled from the host galaxy. This 
would then be compatible with the overwhelming astronomi- 
cal evidence that massive black holes reside at the centers of 
most galaxies. 

Our parameter space is therefore 2-dimensional, 
parametrized by the projections oi, 0.2 of the dimen- 
sionless spins at = Si/M? of the individual black holes 
on to direction of the angular momentum (chosen as the 
z-axis). As a result, spins that are aligned with the orbital 
angular momentum are characterized by positive values of 01, 
a,2, while anti-aligned spins have negative values. Previous 
studies of this parameter space Il2l426ll . have considered the 
recoil velocity and final spin of the merger remnant, and have 
constructed phenomenological formulas for these quantities 
given the initial spins a\ and a2 of the binary. 

In this work, we move our focus to the detectability of 
a given set of binaries in the parameter sub-space of (anti-) 
aligned spins, i.e., for each of these binaries and across a set of 
different masses we calculate the signal-to-noise ratio (SNR) 
for the LIGO JH |H, enhanced LIGOJeLIGO) (23], ad- 
vanced LIGO (AdLIGO) JIJMl Virgo dD], advanced Virgo 
(AdVirgo) JH, and LISApfBl detectors . 

In this way we attempt to address the following questions: 

( i) Which among the aligned-spin configurations is the 
"loudest" and which one is the "quietest"? 

( ii) How large is the difference in signal-to-noise ratio 
between the loudest and the quietest? 
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(Hi) How do these considerations depend on the detec- 
tor used, the mass of the binary, and the number of har- 
monics? 

(iv) Are there configurations whose waveforms are diffi- 
cult to distinguish and are hence degenerate in the space 
of templates? 

Overall, and as expected, we find that equal-spinning, max- 
imally anti-aligned binaries generally produce the lowest SNR 
while equal-spinning, maximally aligned binaries produce the 
highest SNR. For any mass, the SNR can be well described 
with a low-order polynomial of the initial spins p = p(ai, 02) 
and generally it increases with the total dimensionless spin 
along the angular momentum direction, a = |(ai + a 2 ) • L. 
The possibility of describing the whole behaviour of the wave- 
forms from equal-mass, aligned/antialigned binaries in terms 
of a single scalar quantity, namely a, provides a certain 
amount of optimism that also more complex spin configura- 
tions can, ultimately, be described in terms of a few parame- 
ters only. 

We also analyze the impact that modes of the gravitational- 
wave field of order larger than I = 2 but smaller than I = 5 
have on the maximum SNR and show that for low masses 
M £ [20, 100] they contribute, say for the LIGO detector, 
ss 2.5%, whereas for intermediate masses M > 100 Af© they 
contribute w 8% l . In addition, we determine the ratio be- 
tween maximum and averaged SNR for I > 2 which is known 
to be y/E when considering only the I = 2, m = 2 mode. We 
also calculate the mismatch between the waveforms from dif- 
ferent binaries across our spin-diagram and find that binaries 
along the diagonal a\ = —02 cannot be distinguished within 
our given numerical accuracy, whereas configurations along 
the diagonal ai = a-i are clearly different (cf. Fig|7]and[8] as 
well as Table |IVb. Finally, we derive a simple expression for 
the energy radiated in gravitational waves and find that this is 
bounded between ~ 3.6% and ~ 10% for maximally spinning 
binaries with spins anti-aligned or aligned with the orbital an- 
gular momentum, respectively. 

The plan of the paper is as follows: in Sect. H2 we recall 
very briefly the numerical set up and illustrate the properties 
of the initial data used in the simulations. Sect.|III]is dedicated 
to the discussion of the gravitational-wave observables used 
for the subsequent analysis, while Sect. lIVI presents the results 
in terms of the SNR and how this is influenced by higher-order 
modes. This Section also contains a discussion of the match 
between the waveforms from different binaries and an assess- 
ment of the accuracy of our results. Sect.[V] on the other hand, 
provides a brief discussion of the analytic expressions we have 
found representing either the SNR or the energy radiated in 
gravitational waves. Finally, conclusions are summarized in 
Sect ED 



Note that for some specific angles at which the SNR is not maximum, the 
contribution of the higher modes can be much more significant 



II. NUMERICAL SETUP AND INITIAL DATA 



The numerical simulations have been carried out using the 
CCATIE code, a three-dimensional finite-differencing code 
solving a conformal-traceless "3+1" BSSNOK formulation of 
the Einstein equations lfl5ll using the Cactus Computational 
Toolkit 13511 and the Carpet Q6ll adaptive mesh-refinement 
driver. The main features of the code have been presented in 
several papers, and recently reviewed in Pollney et al. lfl5ll . 
The code implements the "moving-punctures" technique to 
represent dynamical black holes following J2, H3l (see also 
JUllHtl), which has proven to be a robust way to evolve black- 
hole spacetimes. 

For compactness we will not report here the details of the 
formulation of the Einstein equations solved or the form of 
the gauge conditions adopted. All of these aspects are dis- 
cussed in great detail in iflal . to which we refer the interested 
reader. More specific to these simulations, however, is the nu- 
merical grid setup. In the results presented below we have 
used 9 levels of mesh refinement with a fine-grid resolution 
of Ax/M = 0.02 and fourth-order finite differencing. The 
wave-zone grid has a resolution of Ax/M = 0.128 and ex- 
tends from r = 24 M to r = 180 M, in which our wave 
extraction is carried out. The outer (coarsest) grid extends to 
a spatial position which is 819.2 M in each coordinate direc- 
tion. Furthermore, because the black holes spins are all di- 
rected along the z-axis of our Cartesian grids, it is possible to 
use a reflection symmetry condition across the z = plane. 

The initial data are constructed applying the "puncture" 
method ll40l - l43ll as described in PPill . We have considered four 
different sequences labelled as "r", "s," "t", and "u" along 
straight lines in the (a±, 02) parameter space, also referred to 
as the "spin diagram" (cf. Table U for details). As shown in 
Fig. Q] these sequences allow us to cover the most important 
portions of the space of parameters which, we recall, is sym- 
metric with respect to the a\ = (12 diagonal. 

We note that similar sequences have also been considered 
in but have here been recalculated both using a 

higher resolution and with improved initial orbital parame- 
ters. More specifically, we use post-Newtonian (PN) evolu- 
tions following the scheme outlined in |45], which provides 
a straightforward prescription for initial-data parameters with 
small initial eccentricity, and which can be interpreted as part 
of the process of matching our numerical calculations to the 
inspiral described by the PN approximations. The free pa- 
rameters to be chosen for the puncture initial data are there- 
fore: the puncture coordinate locations Cj, the puncture bare 
mass parameters mi, the linear momenta p,, and the individ- 
ual spins Si. The initial parameters for all of the binaries 
considered are collected in the left part of Table Q] The initial 
separations are fixed at D = 8 M, where M is the total ini- 
tial black hole mass, chosen as M = 1 (note that the initial 
ADM mass of the spacetime is not exactly 1 due to the bind- 
ing energy of the black holes), while the individual asymptotic 
initial black hole masses are therefore M{ = 1/2. The only 
exception is for the binary s_g, for which D = 10 M. 
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TABLE I: Binary sequences for which numerical simulations have been carried out, with various columns referring to the puncture initial 
location ±x/M, the mass parameters rrii/M, the dimensionless spins 04, and the normalized ADM mass A/ ADM = M ADM / M measured 
at infinity. Finally, the last four columns contain the numerical values of the energy radiated during the simulation using the two methods 
described in the text and the corresponding errors between them, as well as the error to the fitted values. 
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FIG. 1: Schematic representation in the (01, 02) plane, also referred 
to as the "spin diagram", of the initial data collected in Table|I] These 
sequences cover most important portions of the space of parameters 
which is symmetric with respect to the cii — (12 diagonal. 



m. GRAVITATIONAL-WAVE OBSERVABLES 

In this Section we discuss the gravitational-wave observ- 
ables that have been studied from the sample reported in Ta- 



ble U and how these have been used to compute the radiated 
energy, the SNR, the horizon distances and the event rates. 

A. Numerical-Relativity (NR) waveforms 

Although the C CAT IE code computes the gravitational 
waveforms either via the Newman-Penrose curvature scalar 
^4 or via gauge-invariant metric perturbations on a Schwarz- 
schild background, the analysis carried hereafter will be made 
in terms of the latter. While the two prescriptions yield, in 
fact, estimates which are in very good agreement with each 
other and with differences below 2% (see discussion in lfl5lo . 
we have found that the results obtained using gauge-invariant 
quantities have a smaller numerical error, and are thus prefer- 
able. 

More specifically, we compute the gravitational-wave am- 
plitudes h\ m and hf m in terms of the even and odd master 
functions Q\ m and Qf m via the relations ll46ll 

him(t) = hj m (t) - ihf m (t) = Qt m (t) - i f dt>Q* m (t>) , 

(1) 

where the gauge-invariant perturbations are typically ex- 
tracted at a radius of r E = 160A/ (see Sec. |IV D| for a dis- 
cussion of the accuracy of our measurements and ref. lfl5ll for 
a comparison among different extraction radii). 

As mentioned before, all our binaries [but s_g] have ini- 
tial separations of D = 8.0M [D = lO.OAf], which, in 
the parameter space that we have considered, leads to a max- 
imum initial frequency of the numerical waveforms, that is 
^ini = 0.084/M. Depending therefore on the mass M, such 
an initial frequency can be greater than the lower cut-off fre- 
quency of the detector w co for a given source at an arbitrary 
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TABLE II: Initial instantaneous frequencies MtJini and associated 
minimum masses Mmin of the NR waveforms for the different mod- 
els and for each detector according to the corresponding lower cut-off 
frequency (i.e., at 30 Hz for Virgo, at 40 Hz for eLIGO, at 10 Hz for 
AdLIGO/AdVirgo, and at 10~ 4 Hz for LISA). All the values for the 
masses are in units of solar masses. 
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distance. Because we expect that for most masses w co will be 
smaller than cjj n i, we need to provide additional information 
about the gravitational-wave signal in the frequency range be- 
tween uj co and cji n i. This can be accomplished by "gluing" the 
NR waveform with a PN part as discussed in the next Section. 

The values of the initial frequencies and of the associated 
minimum masses M m ; n for each of the detectors considered 
are reported in Table Hfl 

B. Matching PN and NR waveform amplitudes 

The existence of a cut-off mass set by the initial frequency 
of the NR simulations would clearly restrict the validity of our 
considerations to large masses only. To counter this and thus 
include also binaries with smaller masses, we account for the 
early inspiral phase by describing it via PN approximations. 
To produce the PN waveforms, and the PN energy that we are 
using directly in Sec. IV Bl we have used the spinning Tay- 
lorTl approximant used in Hannam et al . 14711. and which is 
based on the PN expressions described in I48U 5 5I1 . The choice 
of TaylorTl is motivated by that fact, that in 114711 it is found to 
be more robust in the spinning case than the TaylorT4 approx- 
imant, which was previously found to yiel d excellent results 
in the nonspinning case IdoI (see e.g., IdoI for a comparison 
of different techniques to obtain the gravitational-wave phase 
information for quasi-circular inspiral). These waveforms are 
3.5 PN accurate in the nonspinning phase, and 2.5 PN accu- 
rate in the spin-dependent terms entering the phasing. The 
gravitational-wave amplitudes, on the other hand, have been 
computed according to ref. l57ll (see also JUt]) to the high- 



est PN order that is currently known for each of the spherical 
harmonic modes that we use. 

A phase-coherent construction of hybrid PN-NR wave- 
forms is rather delicate, and has not yet been achieved for the 
higher spherical harmonic modes we use here (see 0, H for 
some recent work in the case of nonspinning binaries). How- 
ever, for the present purpose of computing the SNR and the ra- 
diated energies, such a construction in the time domain is not 
necessary and all of the relevant work can be done much more 
simply in the frequency domain. In practice, we Fourier trans- 
form the PN and NR waveforms and "glue" them together 
at a suitable "glueing" frequency w g i ue . Since the SNR de- 
pends only on the amplitude of the waveform, [cf. eq. ©], it 
is not necessary to match the PN-waveform in the phase. This 
greatly simplifies the process of waveform matching and ba- 
sically reduces to a simple check of the amplitude matching 
to address the error of the mismatch. Indeed, we have found 
that without any parameter adjustment, the PN-waveform am- 
plitudes match rather well with the inspiral part of the NR- 
waveforms, and result in an error which is usually w 1.5% 
and in the worst case w 4.0% for the binary configuration 
to- The only care which is important to pay in the time- 
domain analysis, and in order to limit the noise artifacts in 
the Fourier-transformed amplitudes, is the use of a window- 
ing function (e.g., a hyperbolic tangent) to smoothly blend the 
waveform to zero before the initial burst of spurious radiation 
and after the ringdown, in order to limit spurious oscillations 
in the Fourier-transformed waveform. A representative ex- 
ample is shown in Fig. [2] where we report the noise strain 
for the Virgo and Advanced LIGO detectors, together with 
the Fourier-transformed amplitude of the PN and NR wave- 
form for the maximally spinning model s§. The waveform 
is assumed to be observed at 6 = 0,0 = for a total mass 
M = 200 M Q and from a distance d = 100 Mpc. The glue- 
ing frequency in this case is at / g i U e = w g i ue /(27r) = 27.14 
Hz. 

Since each £, m mode of the gravitational-wave field will 
have a different initial frequency, we need to make sure that 
they are all properly taken into account when determining the 
glueing frequency, so that, at least in principle 

w g i U e > max(u in i){ m . (2) 

In practice, the initial frequency of our highest mode, £ = 
4, m = 4, has an initial frequency (wmi)44 = 2(wi n i)22 ■ As 
a result, we select the glueing frequency according to the bi- 
nary configuration with the largest initial frequency, i.e., the 
binary to, and take uj g \ uc = 2(w ini ) 2 2 = 0.168/M. We also 
measure how sensitive this choice is, by considering how the 
results are affected when choosing instead cj g i uc ± Acj, with 
Aoj <C Wgi uc . More specifically, for Aw = 0.01 /M we find a 
maximal difference in the computed SNR of ~ 2.0% over 
all configurations and all masses. Note that such a differ- 
ence affects equally the maximum and averaged SNRs (see 
Sect. IIIIDI for a discussion on these two different measures of 
the SNR). Furthermore, a change of Au> in oj g \ uc affects only 
marginally the relative difference between SNRs computed by 
including modes up to I = 2 and I = 4, and also in this case 
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FIG. 2: Noise strain for the Advanced LIGO and Virgo detectors 
and the Fourier-transformed amplitude of the PN and NR waveform 
at 9 = 0,(f> = for a total mass M = 200 M Q at a distance 
d — 100 Mpc for the maximally spinning model sg. The glueing 
frequency is at / g i uc = 27.14 Hz. 



the differences are ~ 2.0%. Overall, therefore, the uncertain- 
ties introduced by the choice of uj g i uc are much smaller than 
the typical error at which we report the SNRs. 



C. Radiated Energy 

Since the total energy must be conserved, we can use the ra- 
diated energy as an important tool to verify the accuracy of the 
gravitational-wave amplitude and thus the overall precision of 
our calculations. More specifically, because it is straightfor- 
ward to determine the initial and the final total mass, it is also 
straightforward to compare the difference in the two with the 
radiated energy. In practice, we compute the initial mass of 
the system as Mini = ^ A dm' while the final mass of the 
merger remnant Mg n is deduced from the properties of the ap- 
parent horizon within the isolated-horizon formalism as first 
discussed in ll59ll and then extensively investigated in l60ll . 
The radiated energy is then simply given by the difference 
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(3) 



and should be equal to the energy that has been radiated 
through gravitational waves during the simulation [46] 
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Overall, we have found that for all binaries the difference be- 



tween E Ta £ and E 



QX. 

ad 



is between — 0.5% and — 4.0% and 



a detailed comparison of the numerical values is reported in 
Table |U In Sect. IV B| we will discuss an analytic fit to the 
computed data that provides a simple-to-use measure of the 
amount of mass radiated during the inspiral, merger and ring- 
down as a function of the initial spins. 



D. SNR, Horizon Distances and Event Rates 



Following ref. 161 
filtering searches as 



we define the SNR, p, for matched- 



P = 77 



matched 



Hf)f 
s h {f) 



df. 



(5) 



where h(f) is the Fourier transform of the time domain 
gravitational-wave signal h(t), defined in the continuum as 



h(f)= / h^e-^dt, 



(6) 



and Sh{f) is the noise power spectral density for a given de- 
tector. Hereafter we will consider the S n (f) for the ground- 
based detectors LIGO, enhanced LIGO, advanced LIGO and 
Virgo, as well as the space-bound LISA interferometer. [The 
associated noise power spectral densities are reported in Ap- 
pendixlAll 

Note that since the SNR © depends on the angle from 
the source to the detector, it is useful to introduce the angle- 
averaged SNR (p 2 ), which can be computed straightforwardly 
after decomposing the gravitational-wave signal in terms of 
spherical harmonic modes. More specifically, using the or- 
thonormality of the spin-weighted spherical harmonic basis 
sYtrn, the "angle-averaged" SNR 



Pavg = (p 2 ) 



df 



s h {f) 



(7) 

can be written as a simple sum of integrals of the absolute 
squares of the Fourier-transformed modes hg m (f) 



(')ll 



df 



S h (f) 



(8) 



and hence it can be evaluated straightforwardly. For each bi- 
nary, distance and mass, we have calculated both the "maxi- 
mum" SNR p max for an optimally oriented detector, i.e., the 
SNR for a detector oriented such that it measures only the 
+ polarization of the gravitational-wave signal, and the aver- 
aged SNR. Here the mass is always meant to be the redshifted 
total mass, i.e., (1 + z)M sourcc , where z is the redshift and 
M sourcc is the mass at the source. For sources at small dis- 
tances, i.e., less than 100 Mpc, then z < 0.024 and hence 
M ~ M source to within a few percent. Identical results would 
have been obtained if we had considered the x polarization. 

It is worth noting that if the gravitational-wave signal is 
modeled simply through the dominant I = 2 = m mode (or 
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in our case via a superposition I = 2 = ±m) 2 , the maximum 
SNR can be deduced from the average SNR after exploiting 
the properties of the spin-weighted spherical harmonic -2Y22 
and _2 Y2-2 > namely, 

Pmax = ^5^ vg (* = 2,m = 2) (9) 

= ^p2 vg (£ = 2,m-±2). (10) 

However, such a relation is no longer true when including 
modes with £ > 2, and the relation between the maximum 
and the averaged value of the SNR can only be determined 
numerically. 

When computing the SNR, a reference distance needs to be 
fixed and we have set such a distance to be d p = 100 Mpc. 
The results of the SNR at d p across the spin diagram can then 
be recast in terms of an "horizon distance", namely the dis- 
tance at which a given binary system with redshifted mass M 
has an SNR equal to a threshold for detectability, which we 
chose to be p = 8, as customary for ground-based detectors. 
The horizon distance is then simply defined as 

d H = d p {^^M VC . (11) 

The quantity dn is clearly equivalent to the SNR but has the 
advantage to provide, at least for detectors not operating at 
large SNRs, a simple estimate of the increase in the relative 
event rate R as 

R~(j^—)\ (12) 

where dft a= _i is the horizon distance of the configuration 
with lowest SNR, i.e., which belongs to the extrapolated case 
a = — 1. Although simple, this formula requires a caveat. Ex- 
pression dTZb is valid as an equality only for small horizon dis- 
tances, namely those for which the redshift is negligible. This 
is because at large redshifts the observed masses would differ 
considerably from the masses at the source. In other words, 
at large redshifts the horizon distances would be different not 
only because of the spin, but also because the masses at the 
sources would be intrinsically different. This clearly impacts 
the deduced event rate as defined in (fT2l , which considers only 
the contributions coming from the spin. Hence, for large red- 
shifts the event rate R defined here serves only as a lower limit 
for masses larger than the optimal one and, vice versa, as an 
upper bound for masses smaller than the optimal. 

To fix the ideas, let us consider a concrete example. Let us 
assume that we have calculated the horizon distance for a bi- 
nary with a = — 1 which, as can be deduced from Fig. [4] and 
will be discussed in the next Section, will lead to the small- 
est SNR for a given detector. We also assume that this binary 



2 Note that in our binary configurations due to symmetry, we always have 

hi m = hi— m 



has a mass at the detector which is smaller than the optimal 
one. Let us now consider a binary with the same mass at the 
detector but with a > —1; this binary will clearly lead to 
a larger SNR but because the masses at the detector are the 
same, the mass of the binary with a > — 1 will be (because 
of the redshift) smaller at the source. As a result, its horizon 
distance will be overestimated, and hence the event rate com- 
ing from ([T2l only an upper bound. A similar argument for 
masses larger than the optimal one would instead lead to the 
conclusion that the event rate R is only a lower bound. 

IV. RESULTS 

In what follows we discuss the results obtained in terms of 
the SNR and how this is influenced by higher-order modes. 
We also discuss the match between the waveforms from dif- 
ferent binaries and an assessment of the accuracy of our re- 
sults. 



A. Horizon distances and SNRs 

The results of the analysis discussed above are nicely sum- 
marized in Fig. [3] which shows the averaged and maximum 
horizon distance dn = dn{a,M) for some of the detec- 
tors considered. As mentioned above, the horizon distance 
has been computed at a reference SNR p = 8.0, and is 
parametrized in terms of the total mass of the system (in so- 
lar masses) and of the average dimensionless spin "a" as pro- 
jected along the orbital angular momentum L 

a = ^(oi + a 2) ■ L = -(ai + 02) ■ e z , (13) 

where L = L/\L\, and the orbital plane has been chosen to 
coincide with the (x, y) plane of our Cartesian coordinate sys- 
tem. More specifically, the top left panel of Fig. prefers to the 
LIGO detector, the top right panel to the Virgo detector, while 
the lower left and right panels refer to the advanced versions 
of both detectors, respectively. 

While quite self-explanatory, these panels deserve some 
comments. First, as expected, the maximum SNR is always 
larger than the average one but the difference between the 
two is not constant, changing both with the total dimension- 
less spin a and with the total mass AI. Second, for any fixed 
value of a, the horizon distance (and hence the SNR) grows 
steeply to a maximum mass and then rapidly decreases to very 
small values of ~ 0(1). Clearly, this reflects the existence of 
a sweet-spot in the sensitivity curve of all detectors. Third, 
for any value of a, the maximum horizon distance/SNR also 
marks the "optimal mass" for the binary M opt , namely the 
mass of the binary whose inspiral and merger is optimally 
tuned with the given detector and hence can be seen from 
further away. Note that the differences between the maxi- 
mum and average SNR are largest in the neighborhood of the 
optimal mass. Fourth, the configuration with spins parallel 
and aligned to the orbital angular momentum are generically 
"louder" than those with spins parallel but antialigned with the 
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FIG. 3: Averaged and maximum horizon distance (1h = dfi(a, M) for the LIGO detector (top left panel), for the Virgo detector (top right 
panel), and for the advanced versions of both detectors (bottom left and right panels, respectively). The horizon distance has been computed 
at a reference SNR p = 8.0. 



orbital angular momentum, with the binaries having a = ±1 
being the "loudest" and "quietest", respectively; this is essen- 
tially the answer to question (i) in the Introduction. 3 . Fifth, 
in the cases of the LIGO and advanced Virgo detectors the 
horizon distance is essentially zero at cut-off masses which 
are ~ 900 M© and ~ 3000 Mq, respectively. Sixth, for any 
fixed value of the total mass, the SNR grows with a and, as we 
will discuss later on, this growth is very well described with 
a polynomial of 4th order (cf. discussion in Sect. IV At . This 
is shown more clearly in Fig. |4] which reports the maximum 
SNR p max for the LIGO detector and for a given set of masses 
at a distance d = 100 Mpc. Note that the growth of p max with 
a becomes steeper for masses M > 200 Mq, for which the 
NR-part of the waveform and hence the plunge and ringdown 
phase dominates. In these cases, the SNR is more then dou- 



3 This behaviour can be easily understood in terms of the orbital dynamics: 
the binaries with larger total angular momentum will have a larger number 
of cycles and hence a larger SNR 



bled between a = — 1 and a = +1. Finally, when going 
from the present LIGO/Virgo detectors to their advanced ver- 
sions, the average horizon distances go from ~ 600/800 Mpc 
to ~ 10 4 /1.2 x 10 4 Mpc, thus with an observational volume 
of the Universe that is increased by a factor of ~ 5000/3000, 
respectively. Note that if we assume a Hubble radius of 
~ 4.1 Gpc, both detectors would effectively detect binaries 
within a large range of masses (e.g., 60 < M/Mq < 500 for 
advanced LIGO) across the whole Universe. 

Figure [3] shows similar information but for the planned 
LISA mission. Since the horizon distance can well exceed 
the whole Hubble horizon, the figure reports the averaged and 
maximum SNR p = p(a 1 M) for sources at d = 6.4 Gpc (z = 
1). Many of the considerations made above hold also for the 
LISA detector, and it is interesting to note that for sufficiently 
high and aligned spins (i.e., a > 0.8), the SNR is > O(10) 
already with binaries having masses > few x 10 3 Mq. 

Finally, the most salient information of Figs.[3]and[5]is col- 
lected in Table [III] which reports the properties of the "opti- 
mal" aligned binaries for the different detectors. More specif- 
ically, the Table reports in its different rows the optimal total 
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LISA mission and for sources at d = 6.4 Gpc (z = 1). 



FIG. 4: Maximum SNR p max = p(a, M) for the LIGO detector 
for a given set of masses at a distance d — 100 Mpc. Note that 
the growth of p max with a is very well described with a low-order 
polynomial which is of 4th order for the optimal mass (cf. discussion 
in Sect. IV At . Note also that the dependence on a becomes stronger 
for masses M > 200 Mq, for which the NR-part of the waveform 
and hence the plunge and ringdown phase dominate. In these cases, 
the SNR is more then doubled between a = — 1 and a = +1. 



TABLE III: Properties of the "optimal" aligned binaries for the dif- 
ferent detectors. Shown in the different rows are the optimal total 
aligned spin a, the optimal total mass in solar masses, the optimal 
maximum p max and average p avg SNRs, the optimal horizon dis- 
tance da (expressed in Mpc and where cH^ 1 is the Hubble radius), 
the lower bound for the optimal relative event rate R, and the glueing 
frequency / g i uo for the optimal binary. The masses have been sam- 
pled with an accuracy of 2.5 Mq for the ground-based detectors and 
of 2.5 x 10 4 M Q for LISA. 





LIGO 


eLIGO 


AdLIGO 


Virgo 


AdVirgo 


LISA 


a 


0.8 


0.8 


0.8 


0.8 


0.8 


0.8 


Afopt (Mq) 


197 


180 


290 


395 


390 


5.35 x 10 6 


Pmax 


87 


175 


1667 


118 


1591 


2.91 x 10 6 


Pavg 


52 


104 


991 


70 


944 


1.77 x 10 6 


d H (Mpc) 


1091 


2190 


> cH' 1 


1476 


> cH- 1 


> cH' 1 


R 


18 


17 


16 


16 


17 


26 


/gluo (Hz) 


27.48 


30.51 


18.71 


13.74 


13.91 


1.0 x 10" 3 



B. Influence of higher £-modes 

As discussed in Sect. IIIIDI it is interesting to consider the 
impact that higher-order modes have on the SNR of equal- 
mass aligned binaries and some representative examples of 
this impact is shown in Fig. [6] The left panel of this figure, 
in particular shows the maximum SNR p max as a function of 
the mass for the highly spinning model sg and for the present 
detectors LIGO and Virgo. Different lines refer to the SNRs 
computed using only the I = 2 multipoles (continuous line), 
or up to the I = 4 multipoles (dashed line). Clearly, the contri- 
bution of the higher modes is most important near the optimal 
mass (i.e., M ~ 200 M for LIGO and M ~ 400 M for 
Virgo) but this is also non-negligible for larger masses, where 
it can produce an increase of ~ 8% in SNR in a detector such 
as Virgo. 

The right panel of Fig. [6] on the other hand, shows the 
ratio between maximum and averaged SNR as a function of 
the total projected spin a for a binary of M = 200 Mq 
(5.35 x 10 6 M B ) an d the LIGO (LISA) detector. As men- 
tioned in Sect. IIIIDI this ratio is not expressed by a simple 
algebraic expression [cf. equation (9)], but needs to be deter- 
mined numerically. Interestingly, this ratio is not constant but 
increases by ~ 10% for larger total projected spins, underlin- 
ing the importance of higher-order contributions as the initial 
spins increase. Overall, therefore, Fig [6] provides the answer 
to question f Hi) in the Introduction. 



aligned spin a, the optimal total mass in solar masses, the opti- 
mal maximum p and average p avg SNRs, the optimal horizon 
distance dn (expressed in Mpc and with H^ 1 being the Hub- 
ble radius), the optimal relative event rate R, and the glueing 
frequency / g i ue for the optimal binary. The masses have been 
sampled with an accuracy of 2.5 M Q for the ground-based de- 
tectors and of 2.5 x 10 4 for LISA. 



C. Match between different models 

A quantity providing a wealth of information is the match 
between the amplitudes of the waveforms from two different 
binaries, so as to quantify the differences in the gravitational- 
wave signal relative to some reference models. The match 
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FIG. 6: Left panel: maximum SNR p max as a function of the mass for the highly spinning model sg and for the present detectors LIGO and 
Virgo. Different lines refer to the SNRs computed using only the i = 2 multipoles (continuous line), or up to the I = 4 multipoles (dashed 
line). Right panel: ratio between maximum and averaged SNR p as a function of the spins ai = for M = 200 Mq (M — 3.53 x 10 6 Mq) 
by including modes up to I — 2 and £ = 4 for LIGO (LISA). In contrast to the case 1 = 2, the I = 4-curve is not constant but depends on the 
initial spins oi, «2 



between two waveforms h\{t) and hi{£) (or a template and a 
waveform) can be calculated via the weighted scalar product 
in frequency space between two given waveforms 



(hi\h 2 



45R 



df 



s h (f) 



(14) 



where hi(f) is the power spectral density of hi(t), the as- 
terisk indicates a complex conjugate, and Sh(f) is the noise 
power spectral density of a given detector. The overlap is then 
simply given by the normalized scalar product 



0[h u h 2 ] 



(h!\h 2 



VWl}(>l2|/l2 



(15) 



Two parameters need to be taken into account when com- 
puting the overlap. The first one is the "time of arrival" tA 
corresponding to an offset in the Fourier-transform of the sig- 
nal exp [iui(t — <a)]- The second one is the "initial phase" <& 
of the orbital motion when it enters the detector band. 

For both of these parameters the overlap should be maxi- 
mized. We have considered two possible ways of doing this. 
The first approach involves the best match, which gives an 
upper bound by maximizing over both of the phases of each 
waveform 



A^best = maxmaxmax{C[/ii, /12]} . 

tA *1 *2 



(16) 



The second way, instead, involves the minimax match, and is 
obtained by maximizing over the phase of one waveform but 



minimizing over the phase of the other 
M 



l mimmax 



maxminmax{0[/ii, ^2]} , (17) 

t A $ 2 



and thus represents a "worst-case" scenario since it gives 
lower matches although one is maximizing over the template 
phase. More details on the maximization procedure can be 
found in l62ll63ll . Note that all the matches computed here- 
after refer to the numerical-relativity part of the waveform 
only. 

A sensible way, if not the most sensible way, of evaluating 
expressions ( fT6] i and ( fTTl ) is to use the binary sq, the nonspin- 
ning binary, as a reference and to compute the overlap with the 
binaries at representative locations in the spin diagram, e.g., at 
the corners for so — ss> so — u 8, s o — or along the main 
diagonal, e.g., s_g — sg. In this way we can assess whether 
the waveform produced by a nonspinning binary can be used 
to detect also spinning binaries and how much the overlap is 
decreased in this case. 

This is shown in Fig. [7] which reports the best and minmax 
matches as a function of mass for a waveform containing only 
the I = 2, m = 2 contribution and refers to the LIGO detec- 
tor. Different lines show the match computed between so and 
other representative binaries, and show the remarkable simi- 
larity between the waveforms of binaries having a zero total 
spin. This is shown by the sq — u$ match, which is essen- 
tially very close to 1 for all the masses considered (cf. also 
Table llVt . This result extends to all the other measured quan- 
tities, such as the radiated energy or angular momentum, and 
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FIG. 7: Best and minmax match as a function of mass for a waveform 
containing only the I = 2, m = 2 contribution and referring to the 
LIGO detector. Very similar behaviors can be shown also for the 
other detectors. 



is not particularly surprising. Indeed, it was already discussed 
by 16311 . although the investigation in that case was restricted 
to what is here the w-sequence. In addition, the equivalence 
between nonspinning binaries and binaries with equal and op- 
posite spins has been exploited in the derivation of expressions 
for the final spin presented in a series of works IU6l - ll9ll . The 
results of Fig.[7]and Table |IV| are therefore a simple example, 
although probably not the only possible one, of a well defined 
region of the space of initial configurations (i.e., those of bi- 
naries with equal masses and opposite spins) which can be 
mapped to an almost degenerate region (i.e., essentially to a 
single point) in the space of templates. This is the answer to 
question (iv) in the Introduction and clearly represents a seri- 
ous obstacle towards a proper estimate of physical parameters 
of the binaries that may be removed, at least in part, only if the 
waveform is measured with a sufficiently high SNR. A proper 
discussion of this problem, as well as the determination of 
other degenerate patches in the space of templates, will be the 
subject of future work. 

An equally remarkable result, presented in Fig. UJ is that 
the overlap is also very high between the nonspinning binary 
and the binary with equal and antialigned spins, sq — S-gl 
also in this case, in fact, the best match is .Mbest ^ 0.9 for 
the range of masses that is relevant here. Slightly smaller and 
decreasing with increasing masses are the best matches com- 
puted when comparing the nonspinning binary with the binary 
of parallel and aligned spins, so that .Mbcst ~ 0.8, but only 
for very large masses. The waveforms appear clearly different 
(i.e., with .Mbcst S 0.6) only when comparing the binaries 



FIG. 8: Best match as a function of the total projected spin a for 
a waveform containing only the i = 2, m — 2 contribution. The 
top/lower panels refers to a binary with a total mass (200/400 Mq) 
which are close to the optimal ones for the LIGO/Virgo or advanced 
detectors, respectively. In both panels the dotted line shows the mini- 
mum best match (0.965) needed for a detection. While the data have 
been computed for the LIGO detector, very similar behaviors can be 
shown also for the other detectors. 



along the main diagonal of the spin diagram, for sg — s_§i al- 
though even in this rather extreme case the differences tend 
to become smaller for smaller masses. Overall, this result 
underlines that even simple waveforms, such as those rela- 
tive to nonspinning binaries, will be effective enough to pro- 
vide a detection for most configurations of equal-mass and 
aligned/antialigned binaries. 

A different way to assess "how different" the waveforms 
are across all of the equal-mass aligned/antialigned spins con- 
figurations considered here is nicely summarized in Fig. [8] 
which shows the best match as a function of the total projected 
spin a for waveforms containing only the I = 2, m = 2 con- 
tribution and referring to the LIGO detector. The top panel, 
in particular, refers to a binary with a total mass of 200 Mq 
that is close to the optimal one for the LIGO/Virgo detectors, 
while the bottom panel refers to a binary with mass 400 Mq 
and close to the optimal one for the advanced LIGO/Virgo 
detectors (cf. Table HVT ). Besides the remarkably smooth be- 
haviour of .Mbcst across all the values of a considered, it is 
clear that the waveform from a nonspinning binary can be ex- 
tremely useful across the whole spin diagram and yield very 
large overlaps even for binaries with very high spins. In both 
panels, in fact, the dotted line shows the minimum best match 
(-Mbcst = 0.965) needed for a detection 16411 . This result is 
reassuring in light of the fact that most of the searches in the 
detector data are made using phenomenological waveforms 
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based on nonspinning binaries. 

For completeness, the results presented in Fig. [7] (as well 
as those in Fig. |9]l are also reported in Table [IVj where the 
different columns show .Mbest and .M m inmax and for wave- 
forms computed either using only the £ = 2, m = 2 contri- 
bution (third and fourth columns), only the £ = 3,?n = 2 
contribution (fifth and sixth columns), or all contributions up 
to £ = 4 (last two columns). Interestingly, the matches among 
the high-order modes, e.g., (s )e=3, m =2 ~ (u s )e=3, m =2, is 
systematically higher than those of the lower ones and remains 
true even for higher modes beyond £ = 3, m = 2 which, how- 
ever, we do not report here. This indicates that in order to 
do high-precision parameter estimation by including higher 
modes it is also important that these modes are accurately re- 
solved, so that they can be clearly distinguished from one an- 
other. 

We generally expect the match to degrade when the wave- 
forms are computed by including higher-order modes (e.g., up 
to £ = 4) and that this degradation will become larger with 
increasing inclination 9. The most notable example is for 
the degeneracy along the diagonal a\ = —a 2 , which should 
be broken by the inclusion of higher-order modes (We re- 
call that these configurations lead to different recoil veloci- 
ties lfl6ll which can only be produced by gravitational-wave 
contributions other than the leading order £ = m = 2 mode). 
For this reason we have computed the sky-averaged match of 
waveforms including modes up to I — 4 (i.e., the "complete" 
waveforms) and the corresponding matches are reported in the 
last two columns of Table [IV] Similarly to what was found 
in ||63ll . we measure a marked decreased in the minmax match, 
but a much smaller decrease in the best match (the latter was 
not considered in ll63ll ). Although our resolution should be 
marginally enough for us to detect such a difference in the best 
match, we also believe that a much higher accuracy is required 
to determine this with certainty. Note also that the matches 
with complete waveforms along other directions, e.g., sq — sg 
or so — s_8 do not decrease and this is simply due to the 
very large mismatch we already have with the £ = 2 = m- 
waveforms (in these cases, in fact, the final black holes are 
considerably different and hence the associated ringdowns are 
expected to be different). 

Finally, we note that although Figs|7]and[8]show data com- 
puted for the LIGO detector, very similar behaviors can be 
shown also for the other detectors. 



D. Accuracy of NR waveform amplitudes 

A reasonable concern that can be raised when looking the 
very high matches between the waveforms in the u-sequence 
is that these are simply the result of insufficient resolution. 
In other words, the waveforms may appear similar simply 
because our resolution is not sufficient to pick-up the differ- 
ences. To address this concern we have computed the overlap 
among the waveforms obtained at three different resolutions 
and for a representative binary with nonzero spins, i.e., tq. 
Clearly, a low match in this case would be an indication that 
our results are very sensitive to the numerical resolution and 
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FIG. 9: As in Fig. [8] but now different lines represent the matches 
obtained when comparing the numerical waveforms of the binary ro 
computed at different resolutions. The matches are computed for the 
LIGO detector, but very similar behaviors can be shown also for the 
other detectors. 



hence the conclusions drawn on the degeneracy of the space 
of templates would be incorrect. 

The results of this validation are presented in Fig.|9]and are 
reported in the last eight rows of Table HVl More specifically, 
shown with different lines in Fig. [9] are the matches obtained 
when comparing the numerical waveforms of the binary r 
computed at low resolution (Ax/M = 0.024) and medium 
resolution (Ax/M = 0.020, which is also the standard one), 
as well as at a medium and high resolution (Ax/AI = 0.018). 
The matches are computed considering only the I = 2, m = 2 
mode and for the LIGO detector, but very similar behaviors 
can be shown also for higher modes or for the other detectors. 

Overall, the results reported in Fig.|9]and in Table IIVI show 

that Afbcst,minmax[Aa;i, Ax 2 \ > -Mbcst,mmmax[ft-1, h 2 ], 

i.e., that the differences we measure in the overlaps among 
two different waveforms hi and h 2 are always larger than the 
differences we are able to measure at two different resolu- 
tions Ax 1 and Ax 2 . In other words, the differences in the 
waveforms across the spin diagram are always larger than our 
numerical errors, even along the degenerate u-sequence (of 
course, as we have a convergent numerical code, the match 
between medium and low resolution is worse than the match 
between medium and high resolution). It is also worth men- 
tioning that as long as the dominant I = 2, m = 2 mode is 
considered, the differences in the matches are well within the 
margin of error for numerical-relativity simulations of black 
hole binaries. A recent work has in fact estimated that the 
differences in the waveforms produced by distinct codes is 
^mismatch = 1 - M « 10" 4 for the last « 1000M of the 
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TABLE IV: Best and minmax matches as computed for the LIGO detector for binaries with different spins in the spin diagram. Different 
columns show Mbest and A^minmax for waveforms computed either using only the £ — 2, m = 2 contribution (third and fourth columns), 
only the I = 3, m = 2 contribution (fifth and sixth columns), or the sky-averaged contributions of all modes up to I — 4 (last two columns). 
Finally the last eight rows show the matches at different resolutions (i.e., Ax/M = 0.024, 0.020, 0.018 or low, medium and high, respectively) 
for the binary ro. 
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0.93695 


0.92143 




200 


90746 


90536 


95647 


0.94521 


89646 


88041 




300 


89491 


89197 

'7 ■ 7* • 1 I ' / 1 


95015 


93814 


87303 


84960 




400 


89369 


89065 


94806 


93550 


85492 


82103 


S_g — Sg 


100 


0.78948 


0.78493 


0.87041 


0.85222 


0.78310 


0.74895 




200 


0.63309 


0.62703 


0.90722 


0.88543 


0.63456 


0.59426 




300 


0.56934 


0.56008 


0.90322 


0.88869 


0.56941 


0.52170 




400 


0.54235 


0.53960 


0.91199 


0.89848 


0.55470 


0.49338 


S_g — Ug 


100 


0.94250 


0.94187 


0.96299 


0.94669 


0.93897 


0.89017 




200 


0.91444 


0.91229 


0.96316 


0.93068 


0.90315 


0.85958 




300 


0.90188 


0.89885 


0.95486 


0.91256 


0.87846 


0.83428 




400 


0.89772 


0.89492 


0.95132 


0.90583 


0.85870 


0.80907 


Sg — Us 


100 


0.87127 


0.86817 


0.87656 


0.84229 


0.85866 


0.80969 




200 


0.79750 


0.79477 


0.83582 


0.81476 


0.79074 


0.73526 




300 


0.74063 


0.73884 


0.83897 


0.80378 


0.73616 


0.68774 




400 


0.71798 


0.71343 


0.84955 


0.81925 


0.71203 


0.66611 


ro 


100 


0.99979 


0.99970 


0.99495 


0.98812 


0.99855 


0.99463 


(0.024, 0.020) 


200 


0.99963 


0.99929 


0.99133 


0.97100 


0.99633 


0.98800 




300 


0.99943 


0.99894 


0.98752 


0.95775 


0.99379 


0.98152 




400 


0.99924 


0.99868 


0.98630 


0.95317 


0.99209 


0.97683 


ro 


100 


0.99990 


0.99989 


0.99873 


0.99299 


0.99881 


0.99639 


(0.020,0.018) 


200 


0.99980 


0.99970 


0.99806 


0.98074 


0.99705 


0.98952 




300 


0.99956 


0.99924 


0.99707 


0.97238 


0.99497 


0.98070 




400 


0.99935 


0.99866 


0.99666 


0.97017 


0.99320 


0.97429 



dominant mode of non-spinning equal mass coalescence l65ll . 
Since the next higher mode I = 3, m = 2 starts to suffer from 
numerical noise, it does not yield the same high agreement, 
and the differences between best and minimax match show a 
larger deviation. 



As a final comment on the accuracy of our waveforms, we 
note that the error made by using waveforms extracted at a 
finite radius, and not extrapolated at spatial infinity is well 
within the error budget of our estimates. We have validated 
this by comparing the waveforms extracted at a finite radius 
against the waveforms computed at future null infinity, via a 
newly developed Cauchy-characteristic code ll6ql . In the case 
of the nonspinning configuration so we have found an error in 
the calculated SNR of less than 1.0% (details on this compar- 
ison can be found in Appendix IBT). 



V. FITTING FORMULAS 

In what follows we provide some simple analytic represen- 
tation of most of the results presented in the previous Sections 
and, in particular, we give a brief discussion of fitting expres- 
sions that can be derived to express the SNR for an optimal 
mass and the energy radiated in gravitational waves. 



A. SNR 

As discussed in Sect. HHP I the maximum SNR depends 
on several factors, most notably on the two initial spins, the 
total mass of the system and, although more weakly, on the 
number of multipoles included in the waveforms. The result- 
ing functional dependencies when one degree of freedom is 
suppressed and the SNRs are presented in terms of the total 
projected spin are shown in Figs. [3] Hand are clearly too cum- 



13 



300 




I i i i i I i i i i I i i i i I i i i i 

-1 -0.5 0.5 1 



TABLE V: Fitting coefficients for the maximum SNR computed for 
the optimal mass [cf. eq. fOfll, The different rows refer to the various 
detectors and have been computed including all modes up to £ = 4. 



detector 


ko 


ki 


k 2 


k 3 


&4 


LIGO 


50.76 


27.11 


13.43 


8.58 


4.63 


eLIGO 


102.45 


53.63 


25.33 


17.67 


11.26 


AdLIGO 


1020.42 


492.25 


243.60 


153.84 


46.99 


Virgo 


71.86 


35.23 


17.140 


10.92 


3.789 


AdVirgo 


968.08 


481.52 


236.45 


140.69 


37.91 



a = — 1 and a = 1 will translate into an increase of a factor 
~ 27 in the event rate. It is therefore likely that many of the 
binaries observed will have high spins and aligned with the or- 
bital angular momentum. This will be particularly true in the 
case of LISA if the prediction that the spins of supermassive 
black holes are aligned with the orbital angular momentum 
will hold full . 



FIG. 10: Different symbols show the numerically computed values 
of p m ax(a, M op t) for the different detectors and represent therefore 
the cross section along the optimal mass of Figs.[3]and|5] Note that 
the SNR for the advanced detectors have been divided by 7 to make 
them fit onto the same scale. 

bersome to be described analytically (although still possible). 

However, most of the complex functional dependence can 
still be captured when concentrating on the best case scenario, 
and hence on the SNRs relative to the optimal mass M opt . The 
behaviour of the SNR in this case is shown in Fig.QJJ where 
the different symbols show the numerically computed values 
of /3max(a, M opt ) for the different detectors. Stated differ- 
ently, Fig. [10] represents the cross section along the optimal 
mass of Figs. [3] and [5] (note that the SNR for the advanced de- 
tectors have been divided by 7 to make them fit onto the same 
scale). 

Clearly, the behaviour of the SNR in this case is sufficiently 
simple that it can be represented with a simple quartic poly- 
nomial of the type 

4 

p max (a; i < 4, M = M opt ) = ^ k n a n , (18) 

whose coefficients k n are reported in Table W\ for the five de- 
tectors considered. 

These results address therefore question ( ii) formulated in 
the Introduction. More specifically, when considering the op- 
timal mass, the ratio of the SNRs for maximally antialigned 
spinning binaries to maximally and aligned spinning binaries, 
i.e., p max (a = l)//9 m ax(a = —1) is ~ 3 for both the LIGO 
and Virgo detectors. This ratio is also preserved when con- 
sidering the advanced LIGO and Virgo detectors. Because the 
event rate scales like the cube of the SNR [cf. expressions (0- 
dl2V l. an increase of a factor ~ 3 in the SNR of binaries with 



B. Radiated Energy 

While the SNR is effectively a measure of the amount of 
energy released during the inspiral, it also incorporates infor- 
mation on the properties of the detectors and is not therefore 
an absolute measure of the efficiency of the gravitational-wave 
emission process. This information can have a number of im- 
portant astrophysical applications, and in particular it can be 
used to study the effect the merger has on the dynamics of 
the circumbinary disk accreting onto the binary when this is 
massive (see 16711 for the first suggestion and 16811 for a recent 
nonlinear study). 

In this Section we present a simple formula to compute 
the amount of energy released and express it only in terms of 
the initial spins. Our formula is restricted to aligned binaries 
and is therefore not as generic as the one recently presented 
in 16911 . which however also requires the determination of a 
larger set of coefficients, some of which have uncertainties of 
~ 100%. As we will show below, the two expressions yield 
results in reasonably good agreement, at least in the part of 
the parameter space we investigate. 

In practice, the expression for the radiated energy -E ra d is 
derived by combining a fit to the numerical data for the bina- 
ries at an initial and finite separation D = 8 M 4 (we refer 
to this energy as to Ef^), with the estimate of the energy re- 
leased from the binary when it goes from an infinite separation 
down to D (we refer to this energy as E™), i.e., 

-Brad = E Tad + £7 ra d = ^ADM — Mft n + , (19) 



Note that for the binary so, we use an initial separation of D = 10M. In 
order to obtain the radiated energy obtained during a simulation starting 
from an initial separation of D = 8M, we only need to recalculate the 
initial ADM mass of the spacetime for this initial separation. The final 
mass of the remnant is in fact the same. 
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where Madm is the initial ADM mass as measured at spa- 
tial infinity of the binary with separation D, and -Mfi n the 
Christodoulou mass of the final black hole 5 . For the fit of the 
radiated energy during the numerical evolution, E™, we use 
the same symmetry arguments first made in lfl6ll and then suc- 
cessfully used in fl7l - [l9ll to write a simple expression which 
is a Taylor expansion in terms of the initial spins 



j%ft( g = l, 0l ,o 2 ) 
M 



Po + a 2 ) + p 2 (a 1 + a 2 f 



Fitting then the numerical data we obtain the following values 
for the coefficients 



Po 



P2 



3.606 ±0.0271 
100 

0.489 ±0.0254 

Too 



pi 



1.493 ±0.0260 



100 



(21) 



where the reduced chi-squared is Xred = 0.008, and where 
the largest error is in the 2nd-order coefficient but this is only 
~ 5%. The different coefficients d2"TT i can then be interpreted 
as the nonspinning orbital contribution to the energy loss (p , 
which is the largest and of ~ 3.6%), the spin-orbit contribu- 
tion (pi, which is < 3.0%), and the spin-spin contribution (p 2 , 
which is < 2.0%). The relative error between the numerically 
computed value of E^ and the fitted one is reported in the 
last column of Table U 

The PN expression for the energy radiated by the binary 
when going from an infinite separation down to a finite one 
r = d, depends on the total mass of the binary, the mass ratio 
and the spin components, i.e., -E™ = E™(r, M, is, oi, az), 
which is the generalization to unequal masses of the energy 
expression used in the definition of the TaylorTl approximant 
in ref. Wfa . However, exploiting the fact that for equal-mass 
binaries the PN radiated energy Ef^ follows the same series 
expansion used for E^, we obtain for M — 1 — q 



M 



= E: 



PN 
rad,0 



+EZ A (ai + a 2 ) + 2 ( ai + a 2 f , 



(22) 



where the coefficients for D = 8 M are given by 

rad '° ~~ 524288 ~ 100 ' 
e p N 985 _ 0.0664 

rad ' 1 ~ 1048576Vf ~ 100 ' 
1 0.00305 



^rad^ 



32768 



100 



(23) 



A rapid inspection of the coefficients (1231 is sufficient to ap- 
preciate that the PN orbital contribution is only ~ 33%, the 



5 Note that Madm + ' s effectively the mass of the system when it has 
an infinite separation. This is approximately set to 1 in most simulations 
but with a precision which is smaller than the one needed here. 



one of the strong-field regime, but also that the spin-related 
PN contributions are mostly negligible, being at most of ~ 4% 
as produced in the last orbits. 

We can now combine expressions d20b-d2TT> with expres- 
sions d22l-(l23]l and estimate that for equal-mass binaries with 
aligned spins the energy radiated via gravitational waves from 
infinity is 



£rad(ai, a 2 ) 



M 



po + pi(oi + a 2 ) +p 2 (ai + a 2 ) 2 , (24) 



(20) where 



Po 



4.826 
100 



Pi 



1.559 
100 



P2 



0.485 
100 



(25) 



Of course these numbers are specific to equal-mass binaries 
and refer to a situation in which the match between the PN 
evolution and the one in the strong-field regime is made at 
a specific separation of D = 8 M. However, we expect the 
results to depend only weakly on this matching separation 
(as long as it is within a PN regime) and hence that expres- 
sions d24"i i and d25l l are generically valid at the precision we 
are considering them here, namely ~ 5%. 

Using expression d24b a number of quantitative considera- 
tions are possible. Firstly, the largest energy is clearly emit- 
ted by equal-mass, maximally spinning binaries with spins 
parallel and aligned with the orbital angular momentum at 
is .E ra d(a = 1)/M = 9.9%. Excluding the astrophysically 
unlikely head-on collision of two black holes moving near 
the speed of light (in which case E Tajd < 14 ± 3% j7(ilo . 
these binaries are among the most efficient sources of energy 
in the Universe. Secondly, equal-mass nonspinning binaries 
lose a considerable fraction of their mass via radiation, with 
E Ta d(a = 0)/M = 4.8%, while maximally spinning bina- 
ries with spins parallel and antialigned with the orbital angular 
momentum have E Ta d(a = — 1)/M = 3.7%. 

Note that expression d24T i is not a strictly monotonic func- 
tion of the total spin and has a local minimum at a\ = a 2 = 
—pi/{^p 2 ) — —0.8 rather than at a\ = a 2 = —1, and yields 
E lad (a = -0.79)/M = 3.6% (cf. Fig. [EE]). Although rather 
shallow, we do not expect such a local minimum. We there- 
fore interpret it as an artifact of the numerical error of our 
calculations (the difference between the energy radiated at 
di = a 2 = — 1 and that at a\ = a 2 = —0.8 is ~ 2% and hence 
compatible with our overall error). Such a local minimum can 
be removed by adding higher-order terms in expression (l20b 
(e.g., up to 4th order in ai+a 2 ) but these improvements are so 
small that they do not justify the use of a more cumbersome 
expression. A comparison between the numerical values and 
the fitting expression l24l is shown in Fig.QTl where crosses 
and squares represent the E^ and E Ta ^ respectively, along 
the diagonal of the spin-diagram (i.e., for ai = a 2 ), while the 
continuous line refers to our fitting expression. Note that such 
a line is a 1 -dimensional cut of a 2-dimensional surface and 
hence it is not expected to exactly fit all points. 

As mentioned above, Lousto and collaborators [69] have 
recently proposed a more general formula that should account 
for the radiated energy in all of the relevant space of param- 
eters, namely for binaries with arbitrary mass ratio, spin ori- 
entation and size. Restricting their expression to the specific 
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10 



where now 




-0.5 



0.5 



Energy radiated during the numerical calculation 



+ .E™ (squares) 



FIG. 11 

(crosses), the total radiated energy iJrad = E™^ 
along the diagonal of the spin diagram, i.e., for ai — a-i. Shown as a 
continuous line is the analytic expressions given here (AEI fit), while 
the dashed line is the one suggested in ref. ll69ll (RIT fit). Note that 
the lines represent 1 -dimensional cuts of 2-dimensional surfaces and 
hence are not expected to fit well all points. Finally, indicated with a 
dotted line is the prediction for the radiated energy coming from the 
point-particle approach of Jsoh and refined in fnfl . 



subset of binaries considered here corresponds to setting in 
their expression (2): Eb = Ee = 0, v = 1/4 and q = 1. The 
resulting expression is then 



1 



1 



1 



m = l Elsco + m E2 + m Es 

+ ^-[E s (a 1 +a 2 )+E A (a 1 +a 2 ) 2 
64 

+ E D ( ai -a 2 ) 2 ] , 



(26) 



where the fitting coefficients have been determined to be E 2 = 
0.341 ± 0.014, E 3 = 0.522 ± 0.062, E s = 0.673 ± 0.035, 
E A = -0.014 ± 0.021, E D = -0.26 ± 0.44 JH, and where 



Eisco 



Vs\ 0.103803 
""3" J + 4 



1 



-{CL\ + 2 ) 



(oi -a 2 f. (27) 



48^ 648^ 
After a bit of algebra we can rewrite d27l ) as 

■ = go+9i(ai+a 2 )+g2(ai+a 2 ) 2 +g 3 (ai-a2) 2 



M 



SO = 


1 / 

A i 1 - 

4 V 


\/8 103803 \ 


Ei 


-r + — a — r 

3 4 J 


lb 


qi = 


1 


E s _ 1.352 




192^ 


+ 64 ~ 100 ' 




fa 


E A „ 
64 ~ " 


0.0219 
100 ' 






5 


E D 0.270 





E3 

64 



5.025 
100 



qs 



2592^2 64 



100 



(29) 



(28) 



Comparing d24l-(l25]l with d28])-(l29l) is now straightforward 
and shows that: the reduced expression from 16911 has a sec- 
ond order contribution ~ (a\ — a 2 ) 2 , which is absent in our 
expression. The remaining coefficients are rather similar but 
not identical. This comparison is summarized in Fig. [TT1 
where the dashed line corresponds to the fitting proposed in 
ref. l69ll . Note that the maximum efficiency for maximally 
spinning black holes predicted by expression ([28} is ~ 8%, 
but our estimate is larger and ~ 10%. Not reported in Fig.fTTI 
is the prediction made in ref. lf72tl . which is linear in the total 
spin and very close to that coming from ( |28| >. 

While the two expressions provide very similar estimates 
for —0.5 < a\ = a 2 < 0.4, they also have predictions dif- 
fering by more than ~ 20% for highly spinning binaries. Be- 
cause both expressions come as a result of a number of sim- 
plifications and assumptions, it is not easy to judge which one 
is the most accurate one, if any. It is useful to bear in mind, 
however, that expressions (l24l-([25l have been obtained from 
a "controlled" set of simulations with small truncation errors 
and therefore have coefficients with error-bars of the order of 
5%. Expressions d28l-(|29]i, on the other hand, because com- 
ing from more extended formulas and thus fitting a wider set 
of different simulations across many groups, have error-bars 
that are intrinsically larger, as high as 100%. In view of this, 
and of the fact that the coefficients are constant, the simula- 
tions carried out here could be used for a new estimate of the 
free coefficients E 2 ,E 3 ,Es, and Ea in (|28| i (Note that be- 
cause in the expression for the radiated energy d20t there is no 
need for a contribution proportional to (ai — a 2 ) 2 , it should be 
possible to set q 3 = and obtain a numerical constraint for the 
presently inaccurate coefficient Ed). Finally, indicated with a 
dotted line in Fig.QT]is the prediction for the radiated energy 
coming from the point-particle approach of ll20ll and refined 
in 

Simulations involving aligned binaries with unequal masses 
will help to settle this issue and provide an extension to our 
expression ( l24l . This will be the subject of future work. 



VI. CONCLUSIONS 

We have considered in detail the issue of the detectability 
of binary system of black holes having equal masses and spins 
that are aligned with the orbital angular momentum. Because 
these configurations do not exhibit precession effects, they 
represent a natural ground to start detailed studies of the in- 
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fluence of strong-field spin effects on gravitational wave ob- 
servations of coalescing binaries. Furthermore, such systems 
are far from being unrealistic and may be the preferred end- 
state of the inspiral of generic supermassive binary black-hole 
systems. In view of this, we have computed the inspiral and 
merger of a large set of binary systems of equal-mass black 
holes with spins parallel to the orbital angular momentum but 
otherwise arbitrary. Our attention is particularly focused on 
the gravitational-wave emission so as to provide simple an- 
swers to basic questions such as what are the "loudest" and 
"quietest" configurations and what is the difference in SNR 
between the two. 

Overall we find that the SNR ratio increases with the pro- 
jection of the total black hole spin in the direction of the or- 
bital momentum. In addition, equal-spin binaries with max- 
imum spin aligned with the orbital angular momentum are 
more than "three times as loud" as the corresponding bina- 
ries with anti-aligned spins, thus corresponding to event rates 
up to 30 times larger. On average these considerations are 
only weakly dependent on the detectors, or on the number of 
harmonics considered in constructing the signal. 

We have also investigated whether these binaries can lead 
to a degenerate patch in the space of templates. We do this by 
computing the mismatch between the different spinning con- 
figurations. Within our numerical accuracy we have found 
that binaries with opposite spins Si = — S2 cannot be distin- 
guished, whereas binaries with spin S± = S2 have clearly 
distinct gravitational-wave emissions. This result, which was 
already discussed in the past l63ll . may represent a serious ob- 
stacle towards a proper estimate of the physical parameters 
of binaries and will probably be removed only if the SNR is 
sufficiently high. 

Finally, we have derived a simple expression for the energy 
radiated in gravitational waves, and find that the binaries al- 



LIGO 



AdLIGO 

S h (f) 
Virgo 



Appendix B: Comparison of waveforms at future null infinity 

A systematic source of error in the results given in this pa- 
per is the finite radius r E = 160 M at which our waveforms 
are computed. In order to determine its influence on the ac- 



ways have efficiencies E va ,&jM > 3.6%. This can become 
as large as E ra( ^/M ~ 10% for maximally spinning binaries 
with spins aligned to the orbital angular momentum. These 
binaries are, therefore, among the most efficient sources of 
energy in the Universe. 
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Appendix A: Sensitivity curves 

For convenience, we report below the sensitivity curves 
used to compute the SNRs that are often difficult to collect 
from the literature. For LISA we we use the same noise curve 
as for the LISA Mock Data Challenge 3 17311 as implemented 
by Trias and Sintes, and made available by the LISA Parame- 
ter Estimation Task Force [74]. The noise curve for advanced 
Virgo can be found in tabulated form in Ref. 1I32I1 . 



curacy of the values reported here, we have exploited the re- 
cent possibility of computing waveforms unambiguously at 
future null infinity J + l66ll . In this approach, which makes 
use of the Cauchy-characteristic extraction technique J75It79Ii . 
the gravitational-wave information ^4 is computed at J + in 
a gauge invariant way and with no causal influence from the 



8h(f) = S {(^£) 56 + 0.16 (£) + 0.52 + 0.32 (^)'j , S Q = 9 x HT 46 , / = 150Hz, 



(X)- 4 - 14 _5(^)% 



Jo t 2U 



- \ 2 



, S = IQ~ 49 , / = 215Hz, 



w) = *{(*SF) 



-4.8 



_6_ 

IT 



So = 10.2 x lO" 40 , f = 500 Hz . 

(Al) 
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outer boundary. 

In practice, we have computed the match between the wave- 
forms extracted at r E and at J + for the nonspinning configu- 
ration so, and found that /Wbcst = 0.999, which is thus within 
the error given by the match between different numerical res- 
olutions (cf. discussion in Sect. IIV and see also Table iTVb. 
Note that the initial separation of the two black holes as re- 
ported in 066I1 . d — 11 M, is larger than the one reported 
here, thus resulting in a much smaller initial frequency u>i n j. 
Nevertheless, we have considered the same glueing frequency 
w giuc = 0.168/M so as to have a fair comparison between the 
two waveforms. 

In addition, we have also compared SNRs obtained in the 
two cases, when the Fourier-transform of h(t) as given in 



terms of '5 4 is easily obtained as 



where ^4 is the Fourier-transform of ^4. For any of the total 
masses considered here and for all of the detectors, we find 
that the differences in the SNRs is less than 1.0%. Overall, 
both results show that the error introduced by the use of a fi- 
nite radius calculation is within our numerical error-bars of 
~ 2.0% and thus does not modify significantly the results ob- 
tained in this work. 
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